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I. INTRODUCTION 



Detection of cosmic rays of energies above 10^° eV has raised yet unsettled ques- 
tions regarding their origin and composition. The first problem is that it is difficult to 
imagine any astrophysical site for the cosmic accelerator (for a review, see Ref. |^). The 
Larmour relation for a particle of charge (E/lO^^eV) = kpc)(|-B|//iGauss), sets 

the scales for the required size, R, and magnetic field strength, \B\, of the accelerator. 
One would expect any sources with sufficient R\B\ to accelerate particles to ultrahigh 
energies to appear quite unusual in other regards. 

A second issue is the composition of the observed cosmic rays. The shower profile 
of the highest energy event is consistent with its identification as a hadron but not 
as a photon [^]. Ultrahigh-energy[| (UHE) events observed in air shower arrays have 
muonic composition indicative of hadrons|l|. The problem is that the propagation of 
hadrons - neutrons, protons, or nuclei - over astrophysical distances is strongly affected 
by the existence of the cosmic background radiation (GBR). Above threshold, cosmic-ray 
nucleons lose energy by photoproduction of pions, A^7 Nir, resulting in the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff in the maximum energy of cosmic-ray nucleons. If the 
primary is a heavy nucleus, then it will be photo-disintegrated by scattering with GBR 
photons. Indeed, even photons of such high energies have a mean-free-path of less than 
10 Mpc due to scattering from GBR and radio photons^]. Thus unless the primary is a 
neutrino, the sources must be nearby (less than about 50 Mpc). This would present a 
severe problem, because unusual sources such as quasars and Seyfert galaxies typically 
are beyond this range. 

However, the primary cannot be a neutrino because the neutrino interaction proba- 

^We use the term ultra-high energy to mean energies beyond the GZK cutoff (discussed below) which 
can be taken to be lO^^-^ eV. 
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bility in the atmosphere is very smaU. This would imply an implausibly large primary 
flux, and worse yet, would imply that the depths of flrst scattering would be uniformly 
distributed in column density, contrary to observation. The suggestion that the neu- 
trino cross section grows to a hadronic size at UHE|Q has recently been shown to be 
inconsistent with unitarity and constraints from lower energy particle physics]^. 

Since UHE cosmic rays should be largely unaffected by intergalactic or galactic mag- 
netic flelds, by measuring the incident direction of the cosmic ray it should be possible 
to trace back and identify the source. Possible candidate sources within 10° of the UHE 
cosmic ray observed by the Fly's Eye [Q] were studied in Ref. 0.0 The quasar 3C 147 and 
the Seyfert galaxy MCG 8-11-11 are attractive candidates. Lying within the la error 
box of the primary's incoming direction, the quasar 3C 147 has a large radio luminosity 
(7.9 X lO^'^ergs"^) and an X-ray luminosity of about the same order of magnitude, in- 
dicative of a large number of strongly accelerated electrons in the region. It also produces 
a large Faraday rotation, with rotation measure RM = —1510 ± 50radm~^, indicative 
of a large magnetic fleld over large distances. It is noteworthy that this source is within 
the error box of a UHE event seen by the Yakutsk detector. However, 3C 147 lies at a 
red shift of about z = 0.545, well beyond z < 0.0125 adopted in Ref. p as the distance 
upper limit for the source of UHE proton primaries. Just outside the 2a error box of the 
primary's incoming direction is the Seyfert galaxy MCG 8-11-11. It is also unusual, with 
large X-ray and low-energy gamma-ray luminosities (4.6 x 10^^ erg s~^ in the 20 — 100 keV 
region and 7 x lO^^ergs"^ in the 0.09 — 3MeV region). At a redshift of 2; = 0.0205, it is 
much closer than 3C 147, but it is still too distant for the flux to be consistent with the 
observed proton flux at lower energies]^. 

Briefly stated, the problem is that there are no known candidate astronomical sources 

^Ten degrees is taken as the extreme possible deflection angle due to magnetic fields for a proton of 
this energy. 
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within the range of protons, neutrons, nuclei, or even photons. Yet there are good 
candidate sources at 100-1000 Mpc. In this paper we propose that the answer to this 
cosmic-ray conundrum may be that UHE cosmic rays are not known particles but a new 
species of particle we denote as the uhecron, U. The meager information we have about 
the cosmic ray events allows us to assemble a profile for the properties of the uhecron: 

1) The uhecron interacts strongly: Although there are only a handful of UHE events, 
the observed shower development and muonic content suggests a strongly interacting 
primary. 

2) The uhecron is stable or very long lived: Clearly if the particle originates from 
cosmological distance, it must be stable, or at least remarkably long lived, with r > 
10^s(mt//3 GeV)(L/l Gpc) where L is the distance to the source. 

3) The uhecron is massive, with mass greater than about 2 GeV: If the cosmic ray 
is massive, the threshold energy for pion production increases, and the energy lost per 
scattering on a GBR photon will decrease. We will go into the details of energy loss 
later in the paper, but this general feature can be understood from simple kinematics. 
In C/7 — > C/tt, the threshold for pion production is Smin = ttiij + rri^ -\- 2'mum.j^. In the 
cosmic-ray frame where the U has energy Eu 3> rriu and the photon has energy ~ 3T 
(where T = 2.4 x 10""^ eV is the temperature of the GBR), s ~ AE^Eu- Thus, 
the threshold for pion production, s > Smin, results in the limit Eu ^ mT^mu / {'^E^) . 
More generally, the threshold for producing a resonance of mass Mji = Mu -|- A is 
Eu = Amu/i^E^). For E-y = 3T, and if the uhecron is the proton, the threshold for 
pion photoproduction is Eu ~ 10^° eV. Of course the actual threshold is more involved 
because there is a distribution in photon energy and scattering angle, but the obvious 
lesson is that if the mass of the primary is increased, the threshold for pion production 
increases, and the corresponding GZK cutoff will increase with the mass of the cosmic 
ray. Furthermore, since the fractional energy loss will be of order m-„/mu, a massive 
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uhecron will lose energy via pion-photoproduction at a slower rate than a lighter particle. 
Another potential bonus if the cosmic ray is not a neutron or a proton is that the cross 
section for f/7 — > f/vr near threshold may not be strongly enhanced by a resonance such 
as A(1232), as when the U is a nucleon. Although there may well be a resonance in 
the Utt channel, it might not have the strength or be as near the pion-photoproduction 
threshold as the A(1232) is in the pion-nucleon channel. 

4) We will assume the uhecron is electrically neutral: Although not as crucial a re- 
quirement as the first three, there are three advantages if the uhecron is neutral. The 
first is that it will not lose energy through e^e~ pair production off the CBR photons. 
Another advantage of a neutral particle is that because it will be unaffected by inter- 
galactic and galactic magnetic fields, its arrival direction on the sky will point back to 
its source. Thirdly, there will be no energy losses due to synchrotron or bremsstrahlung 
radiation. Of course because a neutral particles will not be accelerated by normal elec- 
tromagnetic mechanisms, it is necessary to provide at least a plausibility argument that 
they can be produced near the source. For instance, they may be produced as secondaries 
in collisions induced by high-energy protons. 

In this paper we analyze the possibility that a supersymmetric baryon (wrfs-gluino 
bound state whose mass is 1.9-2.3 GeV - see below) is the uhecron instead of the proton, 
as first proposed in Ref. 0. The has strong interactions, it can be stable, it is 
more massive than the nucleon, and it is neutral with vanishing magnetic moment 0. 
Remarkably, this particle is not experimentally excluded. The light gluino required in 
this scenario would have escaped detection. Experimental limits and signatures are 
discussed in [0 and the reviews of Farrar [§, Q . 

If UHE cosmic rays are S^s, we will show that their range is at least an order of 
magnitude greater than that of a proton, putting MCG 8-11-11 (and possibly even 3C 
147) within range of the Fly's Eye event. 
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While the main thrust of this paper is an investigation into the scenario where the 
is the uhecron, most of our analysis can also be applied to the case where the uhecron 
is much more massive than assumed for the S^. Extensions of the standard model 
often predict new heavy, e.g., multi-TeV, colored particles which in some instances have 
a conserved or almost-conserved quantum number. Bound to light quarks these form 
heavy hadrons, the lightest of which would be stable or quasistable. Such a particle 
would propagate through the CBR without significant energy loss because the threshold 
energy for inelastic collisions is proportional to its mass. Some mechanisms for uhecron 
production discussed below would be applicable for a new very massive hadron. However 
such a particle probably would not be an acceptable candidate for the uhecron because 
its interaction in the atmosphere is quite different from that of nucleons, nuclei, or an 
5°. Although it is strongly interacting, its fractional energy loss per collision in the 
Earth's atmosphere is only of order (IGeV/M), where M is the mass of the heavy 
hadron.0 Thus if the uhecron energy deposition spectrum is indeed typical of a nucleon 
or nucleus, as present evidence suggests, we cannot identify the uhecron with a very 
massive stable hadron. The maximum uhecron mass consistent with observed shower 



properties is presently under investigation 10 



II. PRODUCTION OF UHE S% 



We first address the question of whether there is a plausible scenario to produce 

•^In the infinite momentum frame for the heavy hadron, this is the fractional momentum carried by 
hght partons since they have the same velocity as the heavy parton, but their mass is of order Aqcd- 
It is the momentum of these light partons which is redistributed in a hadronic collision. Of course a 
hard collision with the heavy quark would produce a large fractional energy loss, but the cross section 
for such a collision is small: « a'^/E'^. 
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UHE S^s. This is a tricky question, since there is no clear consensus on the acceleration 
mechanism even if the primary particle is a proton. Here we simply assume that somehow 
UHE protons are produced, and ask if there is some way to turn UHE protons into UHE 
S^s. Our intent is not to establish the viability of any particular mechanism but to see 
that finding a satisfactory mechanism is not dramatically more difficult than it is for 
protons. 

Assuming that there exists an astrophysical accelerator that can accelerate protons 

to energies above 10^^ eV, one can envisage a plausible scenario of production through 

proton collision with hadronic matter surrounding the accelerator. A p-nucleon collision 

will result in the production of -Rp's, the uud-g\umo state whose mass is about 200 MeV 

above the S^. The Rp decays to an and a vr + ,n with the S° receiving a momentum 

fraction of about {mso/mR^y' . From a triple Regge model of the collision, one estimates 

that the distribution of the produced RpS as a function of the outgoing momentum 

fraction x is da jdx ~ (1 — a;)"'^^^"(s')"^~"'^ as x approaches unity. Here, s' = {l — x)s and a 

is the Regge intercept of the SUSY-partner of the Pomeron. Thus, a = ap— 1/2 = e+1/2, 

where e ~ 0.1 is the amount the pomeron trajectory is above 1 at high energies. Hence, 

we parameterize the 5° production cross section in a p-nucleon collision as da/dx = AE^; 

X is the ratio of the energy to the incident energy. Parameterizing the high energy 

proton flux from the cosmic accelerator as dNp/dEp = BE~'^, we have a final fiux 

of dNgo/dE = unLABE^'^^'^ , where nL is the matter column density with which the 

proton interacts to produce an Rp and k is of order 1 (for 7 = 2, k = 0.4). Note that 

the produced S^'s are distributed according to a spectrum that is a bit fiatter than the 

high energy proton spectrum. 

"'The decay Rp S^n was the subject of an experimental search However the sensitivity was 
insufficient in the mass and Hfetime range of interest {m{Rp) = 2.1 — 2.5 GeV, T{Rp) — 2- 10"^" — 2- lO^'^^ 
sec 0) for a signal to have been expected. 
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A disadvantage to this "beam-dump" production mechanism is the suppression 
factor of about AE'^/app, where app is the proton-proton total difFractive cross section. 
This suppression could be of order 10~^ to 10~^ for typical energies. However the pro- 
duced S'°'s enjoy a compensating advantage. The large column densities characteristic 
of most candidate acceleration regions makes it hard to avoid energy degradation of pro- 
tons before they escape. That is, LijipapN + rieCfpe + n^cTp^) may be much greater than 
unity. By contrast, S^'s may escape with little or no energy loss. Their electromagnetic 
interactions are negligible, and analogy with glueball wavefunctions suggests that a so n 
could be as small as IQ^^ap^ 0. Thus the emerging and nucleon fluxes could be of 
the same order of magnitude. This would be necessary for a very distant source such as 
3C 147 to be acceptable, since the required particle flux for the detected flux on Earth 
already pushes its luminosity limit. Assuming that the 3.2 x 10^° eV event of Fly's Eye 
came during its exposure to 3C 147, the resulting time-averaged flux is lleVcm~^s~^, 
which is greater than the X-ray luminosity of 3C 147 |Q. 

In connection with the "beam dump" mechanism, we note that it is possible to have 
at the source a nucleon flux signiflcantly greater than the S*^ flux, and yet at Earth 
still have a large enough flux to account for the high energy end of the spectrum 
without being inconsistent with the rest of the observed cosmic ray spectrum. To see 
this is possible, suppose as an illustrative example that the spectrum for energies 
above lO^^eV is a smooth extrapolation of the proton spectrum at energies below the 
GZK cutoff, i.e., if Jp{E) = AE~^ for E < 10^^-^ eV then JsoiE) = AE^^ for E > lO^o 
eV. Denoting the S'°-to-proton suppression factor by rj, the proton flux for E > 10^^ eV 
is then Jp{E) = rj^^AE^^. Protons of energy greater than the GZK cutoff (here taken 
to be 10^^'^eV) will bunch up in the decade in energy below the GZK cutoff |Tl],|12|. 
The total number in the pile-up region will receive a contribution from protons from 
the source above the GZK cutoff as well as those originally in the pile-up region. With 



rj = 10~^, there will be equal contributions from the pile-up protons and the protons 
originally below the source. The statistics of the number of events with energy above 
]^gi8.5gY poor to exclude this scenario; indeed there is some indication of a bump 

in the spectrum in this region |]T|. 

Note that even for a point source as far away as 1200 Mpc (e.g. 3C 147), the 
required flux of high energy protons at the accelerator is not unacceptable. For instance 
extrapolating the spectrum as 7.36 x lO^^i?"^'"^ eV/m^/sr/sec and using our pessimistic 
efficiency for production (factor of 1/100), requires the high energy proton luminosity 
of the source to be ~ 10^'' ergs/sec. This is indeed a high value, but not impossible. 

Another possible mechanism of high energy production is the direct acceleration 
of charged light SUSY hadrons (mass around 2 — 3GeV), such as Rp and Rn, whose 
lifetime is about 2 • 10^^'^ — 2 ■ 10~^^ sec 0. Due to the large time-dilation factor [E/m k, 
10^^), whatever electromagnetic mechanism accelerates the protons may also be able to 
accelerate the high energy SUSY hadrons. Then, one can imagine that the high energy 
tail of the hadronic plasma which gets accelerated by some electromagnetic mechanism 
will consist of a statistical mixture of all light strong-interaction-stable charged hadrons. 
In that case the flux of the resulting will have the same spectrum as the protons, 
differing in magnitude by a factor of order unity, which depends on the amount of 
SUSY hadrons making up the statistical mixture. Conventional shock wave acceleration 
mechanisms probably require a too long time scale for this mechanism to be feasible 
(e.g., Ref. [|14|)- However, some electromagnetic "one push" mechanisms similar to the 



one involving electric flelds around pulsars [|T^ may allow this kind of acceleration if the 



electric fleld can be large enough. It is certainly tantalizing that the time scale of the 
short time structure of pulsars and gamma ray bursts is consistent with the scale implied 
by the time-dilated lifetime of charged i?-baryons. 

A somewhat remote possibility is that there may be gravitational acceleration mech- 
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anisms which would not work for a charged particle (because of radiation energy losses 
and magnetic confinement) but would work for a neutral, zero magnetic moment particle 
such as an S^. For example, if S'^'s exist in the high energy tail of the distribution of 
accreting mass near a black hole (either by being gravitationally pulled in themselves or 
by being produced by a proton collision), they may be able to escape with a large energy. 
A charged particle, on the other hand, will not be able to escape due to radiation losses. 
Unfortunately, this scenario may run into low flux problems due to its reliance on the 
tail of an energy distribution. 

A final possibility is the decay of a long-lived superheavy relics of the big bang, which 
would produce all light particles present in the low-energy world, including the S^. For 
instance if such relics decay via quarks which then fragment, as in models such as Ref. 
||16|| , the S'''/nucleon ratio is probably in the range 10~^ to 10"^ based on a factor of 
about 10 suppression in producing a 4-constituent rather than 3-constitutent object, 
and possibly some additional suppression due to the S^^s higher mass.0 

Of the scenarios considered above, only the last two are conceivably relevant for a 
super-heavy (0.1 - 1000 TeV) uhecron. Although the energy inp-nucleon collisions {^/s = 



^2Epmp ~ lO^TeV for primary proton energy of lO^^eV) is sufficient for superheavy 



particle production, the production cross section is too small for the "beam dump" 

mechanism to be efficient.^ Also, the direct acceleration mechanism is not useful for a 

superheavy uhecron unless it is itself charged or is produced in the decay of a sufficiently 

long-lived charged progenitor. Even if a sufficient density of superheavy hadrons could 

^After our work was completed, Ref. [ p^ appeared with an estimate of the production of gkiino- 
hadrons from the decay of cosmic necklaces. Note that their pessimism regarding the light gluino 
scenario is mostly based on arguments which have been rebutted in the literature (see for example Refs. 



^ The cross section is proportional to the initial parton density at a; ~ Mu / ^ times the parton-level 
cross section, which scales as M7^^. 




El and @). 
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be generated in spite of the small production cross section, the time scale required for 
the early stages of acceleration could be too long since it is proportional to jS"^. This 
leaves the decay of a superheavy relic (either a particle or cosmic defect) as the most 
promising source of uhecrons if their masses are greater than tens of GeV. 



III. PROPAGATION OF UHE COSMIC RAYS 



To calculate the energy loss due to the primary's interaction with the CBR, we follow 



the continuous, mean energy loss approximation used in Refs. [0 and |T^. In this 
approximation we smooth over the discrete nature of the scattering processes, neglecting 
the stochastic nature of the energy loss, to write a continuous differential equation for 
the time evolution of the primary energy of a single particle. The proper interpretation 
of our result is the mean energy of an ensemble of primaries traveling through the CBR. 
We shall now delineate the construction of the differential equation. 

For an ultrahigh energy proton (near 10^'' eV in CBR frameQ), three main mecha- 
nisms contribute to the depletion of the particle's energy: pion-photoproduction, e+e~ 
pair production, and the cosmological redshift of the momentum. Pion-photoproduction 
consists of the reactions pj — > vr^p and vr+n. Pion-photoproduction, which pro- 

ceeds by excitation of a resonance, is the strongest source of energy loss for energies above 
about lO^OeV, while below about lO^^-^eV, e+e" pair production dominates. For the 
scattering processes (pion-photoproduction and e~^e~ pair production), the mean change 
in the proton energy (Ep) per unit time (in the CBR frame) is 

d£^p(scatter) ■<r^ / \ a r-, / n 

— = — 2^ (mean event rate) x AE (1) 

events 



^Let this be the frame in which CBR has a isotropic distribution. 
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where the sum is over distinct scattering events with an energy loss of AE per event. 
The mean event rate is given by 



mean event rate = —-^f{E^)dE^d$, 



(2) 



where 7 = Ep/rrip is necessary to convert from the event rate in the proton frame 
(proton's rest frame), where we perform the calculation, to the CBR frame, da/d^ is 
the differential cross section in the proton frame,® and / is the number of photons per 
energy per volume in the proton frame. To obtain / we start with the isotropic Planck 
distribution and then boost it with the velocity parameter f3 to the proton frame 
1 



n{E, 



75 ' 



(2vr) 



2E'^ 



(3) 



_exp[7E^(l + /3cos^)/T] - 1 
where 6 is the angle that the photon direction makes with respect to the boost direction. 
Integrating Eq. (j^) over the solid angle^ and taking the ultrarelativistic limit, we find 

1 



/ 



27r27 



1 -exp(-E^/27T)_ 
For AE, the energy loss per event in the CBR frame, we can write 



(4) 



AE{cos 6, pr) = ^rup 



1 + 



■ cos6' 



mr. 



\ 



1 + ^ 



(5) 



where Pr-, which may depend on E^ and cos^, is the recoil momentum of the proton and 
9 is the angle between the incoming photon direction and the outgoing proton direction. 
Putting all these together, the energy loss rate due to scattering given by Eq. @) becomes 
(ii?p (scatter) 



dt 



-7 



-1 



dE,f{E,) 



^The differential is dQdrj [Q and 77 are defined below) for the e+e pair production while it is 
dcosO for the pion-photoproduction. 

^The exact angular integration range is unimportant as long as the range encompasses cos0 = — 1 
(where the photon distribution is strongly peaked in the ultrarelativistic limit) since we will be taking 
the ultrarelativistic limit. 
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where only functions yet to be specified are the recoil momentum and the differential 
cross section (for each type of reaction i). 

For the reaction involving the production of a single pion, the recoil momentum of 
the protons in the proton frame can be expressed as 

Pr{E^, cos 9) = 



2q^E^ cos 9 ± {E^ + rrip) ^AE'^m'^ cos^ 9 - 4mlm.p{Ej + nip) + 
2 [{E^ + mpY - E^ cos2 9] 

where = mplmp + E^) — itl'^/2. When the photon energy E^ is approximately at 

threshold energy of rriT^ + m^/2mp and the proton recoils in the direction 9 = 0, the 

recoil momentum is about 171^^. The recoil momentum is a double valued function, where 

the negative branch corresponds to the situation where most of the photon's incoming 

momentum is absorbed by the pion going out in the direction of the incoming photon. 

Thus, since the positive branch will be more effective in retarding the proton (in the 

CBR frame), we will neglect the negative branch to obtain a conservative estimate of the 

"cutoff" distance. It is possible to work out the kinematics for multipion production, 

but for our purpose of making a reasonably conservative estimate, it is adequate to use 

Eq. (0) as the recoil momentum even for multipion production. 

The pion-photoproduction cross section has been estimated by assuming that the 

s-wave contribution dominates, which would certainly be true near the threshold of the 

production. The cross section is taken to be a sum of a Breit-Wigner piece and two 

non-resonant pieces: 

a(pion) = 2ai^ ® ( " '"'^ " ^ J + ^cTjnultipion 



^"^For example, one can easily verify that the maximum proton recoil during one pion production is 
greater than the maximum proton recoil during two pion production. 

12 



Pcm 



Pcm^X 



tot 



o"nonres 



r(A ^ Xp) 
Ttot 



cTnonres = ^„ \Mip^ ^ np)\' 



Plm 2m^rtot 



^J[m\ - (m^ + nipY] [m\ - {nip - m^) 
1 \J[s - {rup + m^)2] [s - (mp - m^Y] 



v 

^multipion = »tanh (^ ^^ ~^"^^^^^ ) ©(i^, - i^^^iti) (8) 

where lOx is defined through Ahlox = / (ii7|A^(A X'p)\^ ^ Ai denotes an invariant 
amphtude, the center of momentum momentum is given as usual by 



f^p + mxY] [s - [nip - nix 



\21 



pL = f ' ' "^'^ (9) 

and c"ixiultipion ^ crude approximation^ for the contribution from the multipion pro- 
duction whose threshold is at E^^^a = 2(771,^ + m^/mp). The a.„ component of the cross 
section is fitQ to the jry nvr^ data of Ref. , while the amplitude a for o"niultipion 
is estimated from the py Xp data for energies E^ > 0.6 GeV. The numerical values of 
the parameters resulting from the fit are (tu^cj^) = 0.086 GeV^, |A^(p7 ttp)\ = 0.018, 
-'^tot ~ 0.111 GeV, niA = 1.23 GeV, and a = 0.2mb. The factor of 2 multiplying cr^ 
accounts for the two reactions prf vr^p and prf — > vr+n, since a neutron behaves, to first 
approximation, just like the proton. For example, the dominant pion-photoproduction 
reactions involving neutrons are w-f —>■ ir^n and n'-f ir^p which have similar cross sec- 
tions as the analogous equations for protons. Thus, we are really estimating the energy 

^^The functional form was chosen to account for the shape of the cross section given in Ref. pO| . 

^^The fit is quahtatively good, but only tolerable quantitatively. The fit to the data in the range 
between 0.212 GeV and 0.4 GeV resulted in a reduced Xie ^ 50 (due to relatively small error bars). 
This is sufficient for our purposes since our results should depend mainly upon the gross features of the 
cross section. 
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loss of a nucleon, and not just a proton. 

Taking the jry e^e~p differential cross section from Ref. [^] (as done in Ref. 
we use0 



I 47^2 2rjQ SQ^rf rj 2if ) Sml'r]'^ 



1 - 



e: 



Amir] 



2 



+w 



,(10) 



1 /2 

where w = [1 — l/{2Qri — Q'^)] . The recoil momentum is contained in Q = Pr/2me, 
and the photon energy is contained in rj = cos6/2me. 

The final ingredient in our energy loss formula is the redshift due to Hubble expansion. 
We assume a matter-dominated, flat FRW universe with no cosmo logical constant. Thus, 
the cosmological scale factor is proportional to t^/^. The energy loss for relativistic 
particles (such as our high energy proton) due to redshift is then given by 



(iEp (redshift) _ 2Ep 



(11) 



dt 3t 

Furthermore, note that the expansion of the universe causes the temperature to vary 
with time as 

Adding Eqs. (|^) and (|TT]), we have the proton energy loss equation 

dEp (ii^p (scatter) (iii^p( redshift) . , 

~dr ~ Jt ^ Jt ' ^ ' 

whose integration from some initial cosmological time ti to the present time to gives the 

present energy of the proton that was injected with energy Ei at time t^. Note that 

we are interested in plotting Epito) as a function of to — ti with to fixed, which is not 

equivalent to fixing ti and varying to because there is no time translational invariance in a 

^■^We ignore that n does not pair produce e~^e~. However, this has consequences only for energies 
below about lO^^-^ eV. 
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FRW universe. Note also that we need to set the Hubble parameter h (where the Hubble 
constant is lOO/ikms"^ Mpc~^) in our calculation because the conversion between time 
and the redshift depends on h. To show the degree of sensitivity of our results to h we 
will calculate the energy loss for /i = 0.5 and h = 0.8. 

Now, suppose the primary cosmic ray is an instead of a proton. The e^e^ pair 
production will be absent (to the level of our approximation) because of the neutrality 
of S^. Furthermore, the mass splitting between and any one of the nearby resonances 
that can be excited in a jS^ interaction is larger than the proton- A mass splitting, leading 
to a further increase in the attenuation length of the primary. Perhaps most importantly, 
the mass of being about two times that of the proton increases the attenuation length 
significantly because of two effects. One obvious effect is seen in Eq. (^, where the 
fractional energy loss per collision to leading approximation is proportional to py/rrtp 
while Pr has a maximum value of about m^. Replacement of irip m^o obviously 
leads to a smaller energy loss per collision. The second effect is seen in Eqs. (^) and 
(^), where for the bulk of the photon energy integration region, a decrease in 7 (in 
the exponent) resulting from an increase in the primary's mass suppresses the photon 
number. In fact, it is easy to show that if we treat the cross section to be a constant, 
the pion-photoproduction contribution to the right hand side of Eq. (^) can be roughly 
approximated as 

dEJir) m!TV , / 3 4 



- 2 exp(-y/2) 1 + - + - (13) 

dt vr^ \ y y J 

where y = mT^mp/{EpT), clearly showing a significant increase in the attenuation length 
as nip is replaced by m^o. 

The relevant resonances for the S'°7 collisions are spin-1 R\ and i?s [0] (whose con- 
stituents are those of the usual A and S baryons, but in a color octet state, coupled to 
a gluino E^). There are two R-baryon flavor octets with J = 1. Neglecting the mixing 
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Fig. 1: The figures show the primary particle's energy as it would be observed on Earth today if it 
were injected with various energies (10^^ eV, 10^^ eV, and 10^° eV) at various redshifts. The distances 
correspond to luminosity distances. The mass of 5*" is 1.9 GeV in the upper plot while it is 2.3 GeV in 
the lower plot. Here, the Hubble constant has been set to 50kms~^ Mpc~^. 
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between the states, the states with quarks contributing spin 3/2 have masses of about 
385 — 460 MeV above that of the and the states with quarks contributing spin 1/2 
have masses of about 815 — 890 MeV above that of the S^. If we require that the photino 
is a significant dark matter component so 1.3 < Mjio/m^ < 1.6 according to Ref. 



and take the mass of to be about 1.6 — 1.8 GeV as expected, then hes in the 
range 0.9 ~ 1.3 GeV. If we assume that is minimally stable, we have m^o ^ rrip + 
resulting in m^o in the range 1.9 to 2.3 GeV. The other resonance parameters are fixed 
at the same values as those for the protons. 

In Fig. 1, we show the proton energy and the 5*° energy today (with h = 0.5) if it had 
been injected at a redshift z (or equivalently from the corresponding distance]^ with an 
energy of 10^^ eV, 10^^ eV, and 10^° eV. To explore the interesting mass range, we have 
set the mass to 1.9 GeV in the upper plot while we have set it to 2.3 GeV in the lower 
plot. For the cosmic rays arriving with lO^'^eV, the distance is increased by more than 
thirty times, while for those arriving with 10^^'^ eV, the distance is increased by about 
fifteen times. In Fig. 2, we recalculate the energies with h = 0.8. 

Using the mean energy approximation, we can also calculate the evolved spectrum 

of the primary spectrum observed on Earth given the initial spectrum at the source 

(where all the particles are injected at one time). With the source at z = 0.54 (the 

source distance for 3C 147) and the initial spectrum having a power law behavior of E~^, 

the evolved spectrum is shown in Fig. 3. We see that even though there is significant 

attenuation for the number at 3 x lO^'' eV for most of the cases shown, when the 

overall cross section (which was originally estimated quite conservatively) is reduced by 

a factor of half, the bump lies very close to the Fly's Eye event. Moreover, taking the 

Fly Eye's event energy to be 2.3 x 10^° eV which is within la error range, we see that 

-"^"^ Marked are the luminosity distances cIl = H^^q^^ [zgo + (^o ^ 1)(-\/2qo-z + 1 — 1)] where the de- 
celeration parameter go is 1/2 in our ilo = 1 imiverse. 
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Fig. 3: An initial injection spectrum having a power law form of E~'^ is evolved through the 

particle's interaction with the CBR during its 1200 Mpc travel to Earth. The masses of the 5*° and its 
associated resonance are shown. The curve labeled reduced a has the same mass parameters as the solid 
curve except with our conservative estimation of the total cross section reduced by a factor of half. 

the can easily account for the Fly Eye's event. For sources such as MCG 8-11-11, 
clearly can account for the observed event without upsetting the proton flux at lower 
energies. 



IV. CONCLUSIONS 



We have considered the suggestion that the very long-lived or stable new hadron called 
S*^, a uds-ghiino bound state predicted in some supersymmetric models, can account for 
the primary cosmic ray particles at energies above the GZK cutoff. We noted ways that 
conventional acceleration mechanisms might result in acceptable fluxes of high energy 
S'^'s. We also found that the -S'^ can propagate at least fifteen to thirty times longer 
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through the CBR than do nucleons, for the same amount of total energy loss. Thus, if 
S'°'s exist and there exists an acceleration mechanism which can generate an adequate 
high-energy spectrum, 5'°'s can serve as messengers of the phenomena which produce 
them, allowing MCG 8-11-11 Seyfert galaxy or 3C 147 quasar to be viable sources for 
these ultrahigh energy cosmic rays. 

Although much of the relevant hadronic physics in the atmospheric shower develop- 
ment will be similar to that for the proton primaries, some subtle signatures of an S'^ 
primary are still expected. Because an is expected to have a cross section on nucle- 
ons or nuclei somewhere between 1/10 and 4/3 of the p-p cross section, the depth of the 
shower maximum may be a bit larger than that due to the proton. Furthermore, because 
it is about twice as massive as the proton, it deposits its energy a bit more slowly than 
a proton, broadening the distribution of the shower. There may be further signatures in 
the shower development associated with the different branching fractions to mesons, but 
we leave that numerical study for the future. 

A prediction of this scenario which can be investigated after a large number of UHE 
events have been accumulated is that UHE cosmic rays primaries point to their sources. 
If there are a hmited number of sources, multiple UHE events should come from the same 
direction. Also, the UHE cosmic-ray spectrum from each source should exhibit a distinct 
energy dependence with a cutoff (larger than the GZK cutoff) at an energy which depends 
on the distance to the source. The systematics of the spectrum in principle could reveal 
information about both masses of supersymmetric particles and the primary spectrum 
of the source accelerator. 

We noted that the mass range for a new hadron which can account for the observed 
properties of UHE cosmic ray events is limited: it must be at least 2 GeV in order to 
evade the GZK bound, yet small enough that the atmospheric shower it produces will 
mimic an ordinary hadronic shower. 
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